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Abstract. I study the evolution of mean field and linear quantum fluctuations 

C^ ■ in a toroidal Bose-Einstein condensate, whose interaction strength is quenched 

C^ ' from a finite (repulsive) value to zero. The azimuthal equal-time density-density 

OX)' correlation function is calculated and shows temporal oscillations with twice 

the (final) excitation frequencies after the transition. These oscillations are a 

pH , direct consequence of positive and negative frequency mixing during non-adiabatic 

j^ , evolution. I will argue that a time-resolved measurement of the equal-time density 

r^ i correlator might be used to calculate the moduli of the Bogoliubov coefficients 

nH' and thus the amount of squeezing imposed on a mode, i.e., the number of atoms 



excited out of the condensate. 
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1. Introduction 

Bose-Einstein condensates in toroidal traps provide an interesting opportunity to 
study superfluidity and other quantum effects in backgrounds with nontrivial topology 
[IlIlIlllElElElIHllITniiniinilTlIIlIinilTSlin]. These geometries can be created 
{inter alia) by shining a blue-detuned laser onto a harmonically trapped condensate 
thus creating an effectively repulsive Gaussian core. The energy minimum is shifted 
away from the centre of the trap to a finite radius rg typically of order 100 /im 
|181 \W[ [2D]. In present experiments [181 \W[ \^U\ . radial trapping is relatively weak 
so that azimuthal and radial degrees of freedom must be considered in general. Using 
different trapping techniques, e.g., higher Laguerre-Gaussian beams as proposed in 
|21) . tighter radial confinement could be achieved restricting motion in that direction 
to the trap ground state and thus making the system effectively one-dimensional. 

In the azimuthal direction, the condensate is usually (almost) homogeneous and 
obeys periodic boundary conditions permitting stationary solutions with persistent 
current, i.e., non-zero phase gradient, which could be excited through stirring with a 
laser [^ ^^ or orbital angular momentum transfer from a Laguerre-Gaussian beam 
[IHl [ini [21] • By introducing an azimuthal position dependent potential, i.e., putting 
some small obstacles into the torus, inhomogeneous flow profiles can be generated. In 
regions with higher potential, the condensate density will be lower and the velocity 
higher possibly leading to a violation of the Landau criterion and thus an instability 
if the local speed of sound is smaller than the flow velocity. For sufficiently high 
potentials, the barrier can be overcome only through tunneling - a realization of the 
boson Josephson junction, where superfluids with different phases are on either side 
of an (impenetrable) tunnel barrier [2l l3l (61 171 (81 (91 [25] 

Local violation of the Landau criterion can also be understood in the context of 
cosmic analogs. As pointed out by Unruh |26) . sound waves in irrotational fluids 
obey exactly the same evolution equations as massless scalar fields in a certain 
space-time metric, where the effective curvature is generated by the fiuid fiow, see 
|27) for review. The transition from sub- to supersonic flow and vice versa can 
be interpreted as sonic horizons. It becomes thus, in principle, possible to study 
some aspects of cosmic quantum effects in the laboratory, e.g., analogs of Hawking 

radiation |lTl[ll[li[17|[151[Mll3ni[3Ill31[31[31[3Sl[Ml[SZl[3Hl[Mllin], or 

the freezing and amplification of (quantum) fluctuations in expanding spacetimes 
[4a[43l[44l[2a[46l[45l[33[4Il[48l[4g[S0l[5ll[521[53|, where the latter could be achieved 
by varying confinement, thus making the condensate expand |461 [47[ l48l [55) [56] , or 
by changing the two-body interaction strength near a Feshbach resonance [SI], i.e., 
reducing the speed of sound I31[li[nni[13[521[151[Sl[Sl]. 

In this Article, I will concentrate on the latter case and consider the evolution of a 
Bose-Einstein condensate in an isotropic toroidal trap when decreasing the (repulsive) 
two-body interactions. Initially, radial confinement shall be relatively weak, i.e., of the 
same order as the chemical potential, such that more than one radial mode is occupied 
and the system is two-dimensional. An adiabatic reduction of the coupling strength 
would lead to a lower chemical potential, so that the condensate would eventually 
become effectively one-dimensional. For rapid variations of the interactions, however, 
the condensate - classical mean field as well as (linear) fluctuations - might not be 
able to follow these changes; the system would not stay in its ground state throughout 
the evolution. The classical order parameter will be homogeneous in the azimuthal 
direction due to isotropy of the trap, such that any mean field excitations can only 



Time-resolved density correlations as probe of squeezing 3 

be in the radial direction, e.g., as radial breathing motions. Hence, any azimuthal 
excitations must originate from the linear (quantum and thermal) fluctuations such 
that the symmetry of the considered trap allows for an unambiguous discrimination 
between (classical) mean field and linear quantum fluctuations. 

The squeezing of quantum and thermal fluctuations during non-adiabatic 
evolution is a very generic phenomenon occurring in many different physical setting 
such that toroidal Bose- Einstein condensates might serve as quantum simulator, e.g., 
for aspects of cosmic inflation \A2\ |44] . During inflation, the size of the universe 
rapidly increased by a huge factor thus stretching the wavelength of any excitation 
mode. At some point, the quantum fluctuations could not follow the rapid expansion 
any more; they froze and got amplified, an imprint of which is observable as small 
anisotropies in the cosmic microwave background - similar to the small density ripples 
in Bose-Einstein condensates. 

Although in both cases quantum fluctuations get squeezed during non-adiabatic 
evolution, there exist also a few differences between (massless) quantum fields in 
expanding space-times and Bose-Einstein condensates with decreasing interaction 
strength: Most notably, the excitation frequencies go to zero in the former case such 
that the quantum fluctuations truly freeze, whereas the excitation energies are always 
non-vanishing in the latter case and are generally determined by the confinement of 
the condensate, e.g., the torus radius. Since, for non-zero quasi-particle frequency, the 
two-point function of a squeezed state oscillates with twice the excitation frequency, 
such oscillations in the azimuthal density correlations would provide a signature of the 
coherent quasi-particle pairs created during the non-adiabatic evolution. 

To exemplify this, let me consider the squeezed harmonic oscillator: Canonical 
position, X = l/\^2uj{a + a^), and momentum, p = —i\Jui jlia — a^), can be expressed 
through raising and lowering operators, a^ and a, whose time-dependences are just 
given by oscillating phase factors, e.g., a(t) ~ e^*"*a(0). If the oscillator is in a 
squeezed state, the vacuum is not defined with respect to a but rather with respect to 
a different operator 6, i.e., 6|0) = 0, which can be linked to the operators a(0) through 
a Bogoliubov transformation, a = ab + f3*b^ . Obviously, the variances of position 
and momentum now become time-dependent, e.g., (x^) ex |ap -I- |/3p -I- 2\a/3\ cos{2ujt). 
Thus, a time-resolved measurement of (x^)(i) would yield the coefficients from which 
the moduli of the Bogoliubov coefficients, \a\ and \/3\, can be easily inferred and thus 
the amount of squeezing, i.e., |/3p, can be quantified. Similarly, the time-dependence 
of the azimuthal equal-time density correlations might serve as probe for squeezing in 
toroidal Bose-Einstein condensates. Note that this is different from the dynamical 
structure function, which is the Fourier transformed of the two-point function at 
different times. 

This Article is organized as follows: In the next section, I review the field 
equations and their linearization, before introducing quasi-particles in Sec.[3l Density- 
density correlations as observable will be discussed there too with particular emphasis 
on the time-dependence of the two-point function due to squeezing of the density 
fluctuations. The employed trap geometry is briefly review in Sec. 3] before I consider 
two explicit examples in Sec. [S] For realistic trap parameters, the time-dependent 
density-density correlations are calculated for a sudden quench as well as a smooth 
tanh-shaped transition of the interaction strength from a finite value to zero. The 
results will be summarized and discussed in Sec. [6l 
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2. Field equations 

Trapped Bose-Einstein condensate can be described by the interacting Schrodinger 
field equation (in units where h = 1) [601I6T] 



ift* 



--— + V(r) + U(tW'^ 

2 771 



*, (1) 



with external potential V{r) and two-body interaction strength U{t). Well below 
transition temperature, most atoms condense in the lowest state, which acquires a 
macroscopic occupation number. It is therefore convenient to split the field operator 
V? into a macroscopic condensate part, which can be treated classically, and small 
quantum fluctuations [62| 

(vp + X + C) ^ , (2) 



'N 



where N ~ A' A counts the total atom number. The atomic operator A commutes 
with linear and higher-order quantum fluctuations, so that particle number is always 
conserved. In expansion ([2]), the order parameter ^ — 0{\fN) describes the condensed 
atoms, X = 0{N^) are linear quantum fluctuations, and ( ^ 0{N'^) higher orders. 

Setting formally U = 0{1/N), I can expand the fleld equation ([IJ in powers of 
N and obtain the Gross-Pitaevskii equation for the order parameter Vt* [SJ 

idt^ = (-Z- + ^(^) + UimA * . (3) 

\ 2m / 

For the linear fluctuations x follows the Bogoliubov-de Gennes equation [M] 

^dtx = (-^ + V{r) + 2U{t)\^A X + Um'x^ , (4) 

which contains a coupling between x ^^nd its adjoint x^ . The equation of motion for 
higher-order fluctuations C, reads 

^^tC = (^"^ + V{r) + 2U{t)\^\'^ ( + U{t)^^C^ 

+ U{t) {2^x^x + **X' + X^X') + 0[U{t)(\ (5) 

with terms quadratic and cubic in x in the second line acting as source for C,. 
These higher orders C, need to be small for the mean-field expansion ^ to be 
valid. Throughout the rest of this Article, I will not regard Eq. ([5]) any further and 
concentrate on Eqs. ^ for the background and ^ for the linear quantum fluctuations. 
Alternatively to the Gross-Pitaevskii equation ([3]) , the evolution of the condensate 
can be described through Bernoulli and continuity equations for mean field density, 
Qq = I^P, and phase, S'o = arg^. The corresponding linear quantum fiuctuations are 

5q^^*X + ^x\ <55=(l/2zeo)(**X-*X^) (6) 

and kinematics of the phase fiuctuations 5S would, in the low energy limit, obey 
the same evolution equations as a scalar field in a certain curved space-time, see, e.g., 

[liiiTiiiiiiiisiiniiMiiMiisiiMiisiisziisiiMiiiiiiaiiaiTiiii^ - 

tools and techniques from general relativity might be applied. However, in view of the 
factor 1/ Qq appearing in the expression for 5S, linear phase fluctuations are only well- 
deflned for sufficiently large background densities qq - which must evidently fail near 
the boundary of a trapped condensate as the smallness of 5S cannot be guaranteed. 
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Nonetheless, it is still possible to introduce self-adjoint operators 

x+ = ^ (x + X^) , *" " Til ^* " ^^^ ^^^ 

similar to those in Eq. © but without prefactors yfgo since these would eventually 



void linearization of phase and density in regions with small (background) density. 
The evolution equations for x± ^-re most conveniently written by grouping the linear 
fluctuation operators into a two- vector (see also |30| ) 

where the coefficients of the 2x2 matrix V 

A{r, t) = -^— + y(r, t) + 2U{t)\^f - C/(i)5R*' , 
Ztn 

V^ 
I3{r,t) = -^— + Vir,t) + 2U{t)\^\^ + U{t)n-^\ 
Ira 

C{r,t) = t/(t)3^2 (9) 

are self-adjoint Hilbert-space operators [acting on L^(R-^)] depending on trap potential 
V{r,t), interaction strength U{t), as well as classical background ^(r,i). In ^, 5R^^ 
and 3^^ are real and imaginary parts of ^^, respectively. 

3. Quasi-particles and observables 

3. 1 . Quasi-particles 

In order to allow for an unambiguous definition of quasi-particles and thus also the 
vacuum state, I will assume that the condensate is initially in a stationary state, i.e., 
that the order parameter performs trivial oscillations, ^(i) — e^*^*^(0) with chemical 
potential jjl. Quasi-particles can then, in principle, be defined as eigenmodes of Eqs. 
([5]), i.e., by diagonalizing P, see also [SS], though one needs to be careful because of the 
explicit appearance of real and imaginary parts of Vf^ in the field equations: Even if the 
order parameter performs only trivial oscillations, \I'(i) = e^'^*^(0), the coefHcients 
of Eq. ([8]) would be time-dependent and an instantaneous diagonalization of P would 
not yield the proper fluctuation eigenmodes. It is therefore necessary to absorb the 
time-dependence of the initial phase in the definition of the quantum fluctuations 
(which is implicitly done in the fluid-dynamic description) before diagonalizing their 
evolution equations, see |30j . (Note, however, that the order parameter might retain 
a space-dependent phase.) 

Following Ref. [3D], the linear field operators x± can be expanded in terms of 
initial quasi-particle solutions (this expansion can also be derived the other way round 
starting from the eigenfunctions of P |65j ) 

X+(r, t) = UA(r, t)ax + u^(r, t)a\ , 

xAr,t)^X>x{r,t)ax + X>l{r,t)a\, (10) 

where each eigenvalue pair ±A of P is summed only once. The entire space-time 
dependence of x± is now contained in the Bogoliubov functions u\{r,t) and '0\{r,t), 
while the mode operators d^ and d^, annihilating or creating an initial particle with 
energy |A|, are time-independent. The initial values of the Bogoliubov functions ua 
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and X)\ are (for stable modes, i.e., purely imaginary eigenvalues A) given by the right 
eigenfunctions of evolution equations (jH]) , after absorbing the initial phase oscillations 
of the background ^ into the linear operators to render the coefficients of ([5]) time- 
independent. 

3.2. Observables 

The calculation of observables is now straight forward. Using Eqs. (O and ([TU|) . density 
and phase fluctuations ([6]) can be expressed in terms of Bogoliubov functions, ua and 
Vx, and mode operators a\, a\ 

Sg = y2(3fi*x++'5*X-) = \/2(uA5R* + OAS*)aA+H.c., 

SS = -^ (5R*X- - 3*X+) = -i- (0a5R* - Ua3*) ax + H.c. . (11) 

V2g V2g 

Their expectation values (Sg) and {6S) must be zero, as they measure only the 
deviation from the mean, go and ^o, respectively. But their correlations are usually 
non-zero. For the density correlations, I obtain 

{Q{r)g{r')) _ ^ ^ ^ uxir)m{r) + x>xir)^^{r)nlir')mir')+x>l{r')^^i'ir') ^^^^ 



{gir)){gir')) |vl/(r)|2 |vI/(r')P 

and for the phase correlations follows 

£:f^^Qfr.'^\ - 1 ^ 1 vx{r)m{r) - UA(r)SvI/(r) t>l{r')m{r') - ul{r')^^{r') 
b(r)b{r)J - 1 + - ^^-yp ^pyp ,[U} 

where I took the expectation values with respect to the vacuum state |0) defined 
through aaIO) = Voa- Phase correlations could be measured in interference 
experiments after splitting the condensate [SSI IHZ]- Since this might be difficult 
to accomplish due to the nontrivial topology of the toroidal condensate, I will 
focus on density correlations as they can be easily obtained from absorption images 

[MllSilMlllSllZDlinilZllEaElISHllST]. 

3.3. Time- dependent correlations 

Although in Eq. (J12p . any time arguments were omitted, for a gapped excitation 
spectrum, the correlation function usually performs temporal oscillations even after 
the external driving ceased. As will be shown in the following, this residual time- 
dependence can be linked to phase coherence of squeezed excitations and could thus 
serve as probe to detect non-adiabatic evolution. For simplicity, I will assume that 
the condensate settles into a stationary state again after some time tend and thus 
facilitates the introduction of final quasi-particles. Analogous to Eq. (|10p . the linear 
field operators 

X+{r,t)=Uir,t)b, + f:{r,t)bi, 

X-{r,t)=Q,{r,t)b^ + g:{r,t)bi, (14) 

can be expanded in terms of final quasi-particle solutions with creation and 
annihilation operators bj^ and b^^. They are related to initial creation and annihilation 
operators through a Bogoliubov transformation, 6^ = Q^kaOa+Z^ka'^A' which implicates 
for the mode functions Ua = jn^nX + f^/^KA and Ox = 0kQ:kA + 9ll3nX- The density 
fluctuations then assume the form 6g = [/iK(r")e^*l"l*Q!„A + h*^{''')e^^'^^* f3Kx]ax + H.c, 
where the function /i„(r) = fK(r)3ff^l/(r)-|-0K('")S^(r) describes the space-dependence 
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of the final density fluctuation modes but also contains constant factors accounting, 
e.g.; for the initial state, or the vacuum amplitude of the density correlations. Due 
to the mixing of positive and negative frequency solutions as well as the mixing of 
different modes during the non-adiabatic evolution, the coefRcient of the annihilation 
operator has many contributions some of which oscillate with positive frequencies |re| 
and some with negative frequencies — |k| . Hence, the equal-time correlator (summation 
over K, (T, and A) 

+ e-^(l«l+l^l)*/i./i.«.A/3:A + e+»(l'^l+l^l)*/i:/j;/3.Aa;A (15) 

generally comprises many oscillating terms with frequencies ± | A | ± | k | . Bearing in mind 
that h^{r) and h„{r') describe the space-dependence of the final density fluctuation 
eigenmodes, any k, a components of (J15p . e.g., hf^h*^aK.\Oi*^x '^'' ^K.^a'^K.xl^aX^ ^'"^^ ^"^ 
obtained by a suitable (spatial and temporal) projection of the (measured) correlation 
function. After some algebra, the moduli of the Bogoliubov coefficients |q;„a| and |/3„a| 
and the constant contributions to ht^ follow. 

However, since the general calculation is rather tedious, I will discuss the simpler 
case with diagonal Bogoliubov coefficients a^A, PkX oc Sk.\ here. Different eigenmodes 
do not couple and the density correlator ()15p assumes the form 

{Sg{r)6g{r')){t) ^ K{r)hl{r') [\a^\^ + |/3„|2 + 2cos(2|K|t)|a,/3„|] , (16) 

where I omitted the phases of the Bogoliubov coefficients a^, /?« and of the spatial 
function /i„ in the oscillating term (i.e., absorbed them in the exponentials). The 
moduli of the Bogoliubov coefficients Uk. and /3„ as well as the modulus of /i„ (i.e., the 
adiabatic contribution) can be easily calculated from mean value, /Ik/i* (|a„p -I- 1/3^1^); 
and amplitude, 2/1^/1* 10^/3^1, of the temporal oscillations together with the unitarity 
relation ja^p - |/3«,|2 = 1. 

Thus a time-resolved measurement of the equal-time correlation function (|12p 
provides a direct means of detecting non-adiabatic evolution of the linear quantum 
fluctuations, i.e., quasi-particles squeezed out of the vacuum can be observed. Any 
information about the initial state, e.g., whether it is thermal or the vacuum, is only 
encoded in the prefactor hf^ir) such that a clear discrimination between amplification 
and initially present particles is possible - a feature the phonon detection scheme 
proposed in [7S] lacks. Note that the correlation function (|15p is at equal times. Its 
Fourier transformed should therefore not be confused with the dynamical structure 
factor S{k,uj) = / (i^Ma;/(27r)^/2+ie*("*-'='-)(g(r, t)^(0, 0)), cf. e.g. [SHIES], which is 
the Fourier transformed of the density correlations at different times. 

4. Toroidal condensate 

Inspired by recent experiments [THl [ini [20] i I will consider a harmonic trap of frequency 
uj with repulsive Gaussian core of strength Vb and width a 

Vir)^ ^^ + l^oCxp|--|. (17) 

For simplicity, I will assume that the condensate is quasi-two-dimensional with radial 
and azimuthal degrees of freedom. Due to the isotropy of the trap, any mean field 
motion can only be excited in the radial direction, whereas quantum fluctuations 
occur at all azimuthal wavenumbers. After integrating out the radial dependence of 
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the density, g{(j),t) — J drr g{(j)^r,t), the azimuthal density correlations assume the 
form 

{Q{^,t)){g{4>',t)) ^ V2^ gg ^ 2^ ^ g2 ' ^ *^ 

where 

27r oo 

/(• -mufi 

d^ drr— ^[uA(r,t)3fi*(r,t) + 0A(r,t)5*(r,i)] (19) 



are the contributions of the initial quasi-particle modes A, see [Appendix A] for more 
details. 

5. Quench 

In the examples, a condensate of 10^ interacting ^^Rb atoms in a quasi-two- 
dimensional trap with frequency uj — 407rHz, cf. Eq. (J17p . will be considered. The 
Gaussian beam shall have a size of ct = 175 /im and intensity Vq — SOOOfta; so that the 
potential minimum of the torus lies at rmin = 277 /xm. This trap geometry is loosely 
inspired by those used in Refs. [HI [53], though it should be noted that I use a much 
higher intensity Vq of the Gaussian beam and also a lower particle number to facilitate 
a quadratic approximation of the radial trapping with effective frequency ojeff, see 
[Appendix A[ for details. Because confinement in the third direction [z) is supposed to 
be sufficiently tight and enters the calculation only via the effective two-dimensional 
coupling strength [/2D = {2V2n/m)as{t)/aj_, I will not give any particular value for 
uj± nor for the initial s-wave scattering length Cg. However, the general features, i.e., 
the oscillations of the azimuthal density correlations, should be retained for three- 
dimensional trap geometries, higher particle numbers, or less intense Gaussian cores 
as well. 

The interaction strength is quenched from Uq ~ 0.025fia;offaeff to Ui = 0. The 
initial chemical potential is //in = 7.29 Woff, so that the lowest 9 radial trap modes 
will have significant occupation numbers and the condensate is thus in a quasi- 
two-dimensional regime. When adiabatically lowering the nonlinear coupling U, the 
chemical potential would decrease and the population of all radial trap modes except 
the lowest would decline - the system would become effectively one-dimensional. For 
rapid changes of U, this, however, is generally not the case any more and (radial) 
breathing oscillations of the mean field will be excited. Several radial trap modes will 
have macroscopic occupation. Evidentially, as the final state is not stationary, the 
definition of a chemical potential is not meaningful any more. 

The results presented in the following have been obtained using a basis expansion 
of the order parameter, XP, and the linear fluctuations, x±, see [Appendix A] and 
keeping the lowest 15 radial basis functions. The initial state was propagated using 
a fourth-order Runge-Kutta algorithm with adaptive step width for the coupled 
evolution equations of mean field ()A.2[) and fluctuations (JA.3P . 

5.1. Sudden quench 

A sudden reduction of the interaction strength U from a flnite value to zero is probably 
the easiest quench dynamics that can be realized in experiments. In Figure (TJ the 
time-dependence of the azimuthal density correlation spectrum is shown for such a 
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Figure 1. Fourier components |5p„p/£iQ of the density correlator as a function 
of azimuthal mode n and time t for a sudden quench of the interaction strength 



from Uq = 0.025 hw^^a^ff to C/i = after fg 



Is. 



0.025 



0.025 r 




Figure 2. Density-density correlations for sudden quench from Uo = 
0.025 hw^ga^g to t/i = at tsweep = 1.5. On the left, the spectrum is plotted 
at different times, whereas the time dependence of the lowest six components 
I'^fiP/fo *° l'5£?6p/£?o i^ shown on the right. 



transition from C/q = 0.025 /k^Jcffaoff to Ui = 0. Relative density correlations on the 
percent level can be observed for low azimuthal excitations. However, the n = 
mode is not squeezed as this would correspond to fluctuations in the (conserved) total 
particle number. The left plot of Figure [5] shows the spectrum at different times, 
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while the right plot displays the time-dependence of the the lowest azimuthal modes. 
The correlation spectrum does not converge to a constant value, but rather exhibits 
temporal oscillations well after the quench even though the interaction strength has 
been tuned to zero. 

These temporal oscillations of correlation functions might seem a bit puzzling at 
first but are actually a direct consequence of the non-adiabatic evolution during the 
interaction quench together with a gapped final excitation spectrum, see also Section 
[321 especially Eqs- (fTS]) and (fT6|l . The oscillation frequencies observed in the plots are 
in good agreement with the expected values: After the interactions have been tuned to 
zero, the azimuthal excitation frequencies are a;„ = n'^ /2mr^^-^^ = n^ x 0.306 Hz above 
the corresponding mean field in each radial mode. For n = 4, 7i = 5, and n = 6, this 
implies period times T4 = 6.42 s, T^ =4.11 s, and Tg = 2.85 s - the values observed in 
Figures [T] and [H For lower excitations, one oscillation period lasts much longer such 
that no full cycles of |5£'nP(i) can be observed during the propagation time. Whereas 
higher excitations undulate very quickly, e.g., n — 20 has repetition time T20 = 0.26 s, 
making it much harder to resolve in Figure [T] due to finite grid size of the plot. Also, 
the amplitude of the oscillations declines with n, as it is related to the Bogoliubov |/3| 
coefficients, which, for a sudden transition, read 

l/3| = I 



^out 



(20) 

'^out 

The azimuthal Bogoliubov coefficients, i.e., summed over all radial excitations, can be 
calculated from minimal and maximal values of the oscillations in Fig. [2] (though it 
might be difficult to extract these values from the plots). From the numerical data, 
I have for the mode n = 4 as minimum |i5p4|^;jj/pQ = 2.44 x 10^^ and as maximum 
l^£'4lmax/£'o = l-"^^ X 10^^, such that the Bogoliubov coefficient of that azimuthal 
mode is |/3„=4p « 18. This relatively large number can be understood from the 
huge energy difference before and after the quench: Initially, the condensate is in an 
interacting regime and the quasi-particles are phonon-like excitations whose frequency 
is dominated by the radial contribution of order of the effective trap frequency, 
oJeS = 281 Hz, cf. [Appendix A| After the quench, however, only the small azimuthal 
C(l Hz) part remains as the condensate becomes non-interacting. 

The oscillations in the azimuthal equal-time density correlations can be 
interpreted as squeezing of quasi-particle (atom) pairs out of the quantum vacuum. 
Due to angular momentum conservation, one of the partners must have azimuthal 
wavenumber +n and one — n, i.e., one moves clockwise and one counter-clockwise. 
If both quasi-particles are at the same place, they will not contribute to the non- 
local correlations, hence the minima of (|16p. As they are moving away from each 
other, their contribution to {g((f>)g{(j}')) grows until it reaches its maximum when both 
quasi-particles are farthest apart. A time-resolved measurement of the density-density 
correlations can be used to trace these coherent pairs and thus infer the amount of 
squeezing incurred during the transition, e.g., the Bogoliubov coefficients, |/3riP, giving 
the number of created quasi-particles in each azimuthal mode (summing over the 
different radial components) |43j . 

5.2. tanh(7i) sweep 

As second example for the dynamics, I will consider a smooth transition (7 > 0) 

Uit) = ^°"t^' - ^°o^' tanh [7(t - Wop)] (21) 
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Figure 3. Time-dependence of ttie correlation spectrum for a smooth tanh- 
shaped transition between Uq = 0.025 TicJoffacff and Ui = 0, cf. Eq. II2H . The 
transition occurs at tsweep = 1 s with sweep rate is equal to the harmonic trap 
frequency 7 = 



"-sweep 

:407rHz. 
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Figure 4. Density-density correlations for a smooth tanh transition from 
Uo = 0.025 tko^ga^ff to C/i = with sweep rate j = uj = 407rHz at tgwoop = Is, 
cf. Eq. 1 I2II I. The left plot shows the spectrum at different times and the right 
shows the time-dependence of the lowest azimuthal correlation. 



from Uo ~ 0.025 hujcBO-cS to Ui = 0, i.e., between the same limiting values as before. 
The sweep rate 7 shall be equal to the trap frequency w, such that only the low- lying 
excitations with energies of about the same order as the trap frequency will evolve 
non-adiabatically. The results are plotted in Figures [3] and |4l Again, correlations on 
the percent level can be observed, though, as expected, smaller than in the previous 
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example. Temporal oscillations of |(5g^j|, characteristic for non-adiabatic evolution, can 
be observed for low excitations only and subside much more rapidly with increasing 
n than for a sudden transition. The spectrum |(5g„p becomes (almost) constant in 
time for higher-energetic modes n w 15. ..20 already because the excitation energies 
are smaller than the inverse sweep rate I/7 and these higher fluctuations can (at least 
partially) adapt to changes of the background. 

6. Discussion 

In summary, I studied the evolution of quasi-two-dimensional Bose-Einstein 
condensates in isotropic toroidal traps when tuning the (repulsive) interaction strength 
from a finite value to zero. Since present experiments [THl[Tni[ini provide only relatively 
weak radial confinement, the initial interaction strength was chosen such that the 
condensate is in an effectively two-dimensional regime with radial and azimuthal 
degrees of freedom. In the third direction (z), confinement shall be sufficiently tight 
such that motion is restricted to the trap ground state and can be integrated out. 

Due to the isotropy of the trap, the classical mean field must be homogeneous 
in the azimuthal direction. Hence, any background motion must occur radially. 
The linear quantum fluctuations, on the other hand, are only subject to (angular) 
momentum conservation and might thus have non-trivial azimuthal dependence, which 
can be treated through a Fourier transformation. For each Fourier mode, two coupled 
one-dimensional field equations (in the radial direction) with gap n^ /2mr^^^ follow, 
which can be easily solved numerically. By integrating out the radial dependence and 
considering only the (relative) azimuthal correlations, a clear signature of the quantum 
fiuctuations can be provided. 

The radial dependence of mean field and linear quantum fluctuations was 
expanded in terms of harmonic oscillator functions. For the initial flnite interaction 
strength, the condensate (mean fleld) is radially distributed over several modes. When 
slowly (adiabatically) reducing the interaction strength U, more and more atoms 
would gather in the trap ground state, while the occupation of all higher modes 
would diminish until the condensate becomes effectively one dimensional. For rapid 
variations of U, however, the condensate might not be able to follow these changes of 
the interaction strength and will therefore not stay in its ground state. A fraction of the 
atoms would remain in higher trap states and the mean field would split up into several 
parts, whose phases oscillate with different frequencies - the condensate would undergo 
classical oscillations in the radial direction. It should be noted that, even though the 
condensate performs breathing motion and is thus split into different parts evolving 
with vastly different frequencies, the mean field is still coherent, i.e., the phases of 
the classical excitations and the non-excited bulk are still uniquely related. Coupling 
to an environment or the back-reaction of quantum fiuctuations might, however, lead 
to loss of phase coherence so that the condensate becomes fragmented into several 
mutually incoherent parts [751 ES] ■ 

For non-vanishing interaction strength after the quench, U ^ 0, and if the 
mean field undergoes radial breathing, the background-dependent terms in the linear 
field equations, C/|\I'(r, i)p and C/^^(r, i), would become periodic in time and could 
amplify quantum fluctuations through resonance. Therefore, I considered only the 
case, where U is tuned to zero and the theory becomes non-interacting, i.e.. the 
linear excitations evolve independently of the condensate background. In view of 
the azimuthal homogeneity of the mean field, azimuthal density correlations - the 
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radial dependence has been integrated out - represent a suitable observable for the 
quantum fluctuations independent of any background motion. Two different time- 
dependences of the interaction strength were considered: firstly, a sudden change from 
a finite C/q to zero and, secondly, a smooth tanh-shaped transition between the same 
values of the interaction strength. In both cases, squeezing of quasi-particles could 
be observed, which manifests as temporal oscillations in the azimuthal correlation 
functions. This can be understood in a simple picture: due to angular momentum 
conservation, quasi-particles are always excited as pairs with azimuthal wavenumbers 
±n. With one quasi-particle (atom) circling the torus clockwise and the other counter- 
clockwise, the oscillations in the correlation function can then be understood through 
co-incidence of the pair. For a sudden transition, all modes evolve non-adiabatically, 
whereas a smooth tanh shape of the interaction strength violates adiabaticity only for 
low excitations, while higher modes stay closer to their ground state. In both cases, 
relative (vacuum) density-density correlations on the percent level could be observed. 

Temporal oscillations in the equal-time correlation function after the quench are 
closely related to finite quasi-particle frequencies and non-adiabatic evolution of the 
fluctuations - initial particle solutions oscillating with positive frequency before the 
quench usually comprise positive and negative frequency parts afterwards, where the 
amount of mixing can be quantified through Bogoliubov coefficients a and /3. For non- 
zero /3, i.e., non-adiabatic evolution, the density correlator will assume oscillations 
with twice the final quasi-particle frequency. From mean value and amplitude of 
the oscillations, the Bogoliubov coefficients as well as the adiabatic correlations can 
be obtained. Thus, a time-resolved measurement of the equal-time density-density 
correlation function permits the determination of the Bogoliubov coefficients |/3p as 
signature of quasi-particle squeezing. 

In contrast to methods aiming at the detection of single phonons, e.g., |75j . the 
temporal oscillations of the density correlation are as signature of squeezing relatively 
robust against changes of the initial quantum state. For instance, a thermal occupation 
number would merely appear as prefactor to the correlations of a mode, whereas the 
Bogoliubov coefficients \(3\ as measure for the amplification of the initial fluctuations 
(and quasi-particles) follow from the relative amplitude compared to the mean of the 
oscillations, which is unaffected by such a prefactor. If the excitation frequencies 
would go to zero, which happens in expanding Bose-Einstein condensates after the 
trap has been turned off |461 \55\ [56l [58l [57] or during cosmic inflation |421 144] , the 
density correlations would converge and a time-resolved measurement would yield no 
additional information. In particular, it would not be possible to distinguish between 
contributions stemming from initial particles and those due to the amplification of 
fluctuations, i.e., squeezing. 

A time-dependent measurement of the equal-time density correlations is also 
fundamentally different from obtaining the dynamical structure factor [60l [76] , which 
is the temporal Fourier transformed of the density correlator at different times. To 
exemplify these differences, let me consider a simple one-mode model g = f{t)a -f 
f*{t)a^ . Assuming that the system is in squeezed state, f{t) = ae^*"* + /?*e*"*, 
the propagator reads {git)git')) = f{t)f*{t') = (jap + a/3e-2''^*')e-''^(*-*') + (|/3|2 + 
a*/3*e^*"* ■jgia;(t-t ) ^^^^ |-jjg dynamical structure factor has two peaks at ±w with time- 
dependent coefficients. It is, in principle, possible to obtain the Bogoliubov coefficients 
and thus the amount of squeezing from the time-dependence of these peaks as well, 
though a time-dependent measurement where both times t and t' are varied would be 
necessary. 
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Appendix A. Basis expansion 

Appendix A.l. General basis 

Method of choice for the numerical propagation of mean field \E' as well as quantum 
fluctuations x± is a basis expansion. To remain general, I will formally consider 
non-orthogonal bases in the following, i.e., tensor notation will be adopted and I will 
distinguish between covariant {/^(r)} and contravariant {f"{r)} bases. They are 
normalized {f°',fi3)L^ = 5°" p with L^ inner product {f,g)L^ ~ J d^rf*{r)g(r). Non- 
orthogonal basis might occur, e.g., when using the eigenfunctions of a non-symmetric 
operator as basis, cf. also 'D |65j . or when picking a particularly suitable numerical 
basis. 

The components of the order parameter and the linear fluctuations 

*(i,r) = Ur)^^{t), x±{r,t) = /a(r)x^(t) (A.l) 

follow from their projection onto the covariant basis ^" = (/",^)^2 and x± = 
(/"iX±)l2. The Gross-Pitaevskii equation ([3]) reads in basis expansion 

i9t*" = /C"/J^^ + C/X"^^^(*'^)**'^*^ (A.2) 

where the entire space-dependence (and possibly the time-dependence of the external 
confinement V) is integrated out and thus contained in the coupling tensors /C"^ = 
J d"r{f°')*{~V^/2m + V)fp and M^'p-ys = J d^r {f"')* f*f^fs. For the linear 
fluctuations, I obtain 



The coupling tensors A°'p = Jd^r{f")*Afp, i3"/j ^ J d^r {f")*Bfp, and C^p = 
/ d^r (/")*C//3 follow by taking the matrix elements of operators ([9]) with respect to 
the basis {/«} and its dual {/"}. Note that, due to the non-linearity of the original 
field equation ([Ij and thus the appearance of the mean field '$ in the linear evolution 
equations ([5]) for x±, the matrix coefficients in (|A.3p usually contain integrals of four 
basis functions similar to M"'p-ys and summation over the components of the order 
parameter ^P. 

From Eqs. (fTO)) and (jA.ip . I can infer the basis expansion of the Bogoliubov 
functions ua and Oa 

UA(r,i) = Ur)ut{t), t,A(r,i) = UrMit) ■ (A.4) 

so that the linear field operators x± can be described by time-independent operators 
a/^, a A creating or annihilating an initial quasi-particle with energy |A| and time- 
dependent coefficients u"(i) and t>A(^)- ^^^ initial values of the Bogoliubov functions 
are given by the (up and down) components of the right eigenvectors of the coupling 
matrix appearing in (|A.3|) . and can be propagated using Eq. (jA.3p . As the evolution 
equation is first order in time, positive and negative frequency solutions appear 
separately in the spectrum and one needs to be careful not to double count the 
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pair, i.e.. (for stable modes) eigenvalues with negative imaginary part correspond 
to annihilation operators a\ and those with positive imaginary part yield creators 
a\. Inserting expansions (|A.1|) and (|A.4|) into ((T^ . observables such as the density 
correlations can be calculated. One should, however, bear in mind that the basis 
functions fa are usually not real and thus complex conjugation of the coefficients u" 
etc. does not yield the coefficients of (u")* etc. with respect to the basis {fa} but 
rather those with respect to {/„}. 

Appendix A. 2. Toroidal condensates 

Toroidal Bose-Einstein condensates can be created using potential p?| . which has its 
minimum at r^^^^ = cr^ \n{2Vo/mLu'^a'^) > for 2Vo > ruLo'^a^ . In order to simplify the 
analysis a bit, this potential can be expanded to second order about this minimum 

V{r) = V^(r„,i„) + ^(r - r^.^f + 0[{r - r^.^f] , (A.5) 

and an effectively harmonic potential with frequency w^g — 2lo'^ ln(2Vo/TO(u;^cr^) = 
2u;^r^ijj/cr^ can be obtained. The first term, V{r^in), yields only a constant energy 
shift and will be omitted in the following. This expansion (jA.5|) of the potential can 
be used if the radial extend Ar of the condensate is much smaller than the torus 
radius rmin, which can be achieved by demanding that the harmonic oscillator length 
Ooff = 1/y'muJcff be much smaller than r^in and only the lowest few radial modes 
being occupied, i.e., a chemical potential of order coes. 

In view of the symmetries of the trap - the potential depends only on r = \r\ and 
is independent of the azimuthal angle - it is advantageous to change to planar 
polar coordinates, r = rer{(j)), and to introduce bases for radial and azimuthal 
dependence. Harmonic oscillator functions ha{r — Tmin) centered at r^in seem to 
be an obvious choice for the radial basis. But one should, however, bear in mind 
that polar coordinates yield an additional factor r as measure in the inner product 
J d^rhahp — j rdrd4>h*ahp. It is therefore more convenient to absorb this factor and 
to use slightly different functions 

/a^r-rmin) = p , (A.6) 

which are approximately orthonormal J„ drr fa{r — rniin)//3('' — ''min) = ^ap- 
Corrections to this orthonormality are exponentially suppressed ex e~'-"+'^-'''°/''"'". 
[Of course, this suppression does not hold for higher modes a > J'min/aoi but the 
occupation numbers in these modes are very small because the condensate is radially 
localized near rmin- Also, expansion (jA.SP of the radial potential ([17]) to second order 
would break down first.] 

In view of the isotropy of the potential (J17p . the azimuthal dependence of mean 
field ^ and linear fluctuations x± is most conveniently described in terms of plane 
waves 

g™ ((/.) = ^e™^ meZ (A.7) 



/27r 

obeying periodic boundary conditions, gm{4>) = 9m{4> + 27r), and normalization 
L d(f>g*^{(j})gn{(t>) = Snm- E)ue to isotropy of the trap (|17p . the classical mean 
field, \I'(r, ^, i) = -^ fa{r)'^°' (t) , must be independent of the angle (j) and the 
azimuthal label on \|>™=''^" can be omitted. The Bogoliubov functions u\{r,t) — 
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9m{<P)fa{r)ul'-°'{t) and Ox{r,t) = gni{(p)fa{ry^'"{t), however, depend on (p and can 
be propagated through Eq. (jA.3|) . Different azimuthal modes decouple (except, of 
course, ±n), because the background does not depend on 4> so that for each angular 
basis function, (?„((/)) ex e*"''^, a set of coupled evolution equations follows. 

Since I consider fluctuations above a dynamical background, I must choose 
observables affected as little as possible by any mean field motion. For any changes 
of the interaction strength, radial breathing oscillations of the order parameter ^ 
are usually excited, whereas the background stays at rest in azimuthal direction. I 
will therefore pick azimuthal density correlations in the following, where the radial 
dependence has been integrated out 

m = Jdrrg^ ^+ g„.meTax + g*^m6eTr&l (A.8) 

The measure r has been included for convenience. The mean azimuthal density, 
N/2tt, is constant due to isotropy of potential (|17p . while the coefficients 5gJ^ = 
(1/v^) I]a(uA'"^*" + 0™"S*") of the linear density fluctuations give the nontrivial 
contributions of the initial eigenmode A to the density fluctuations with azimuthal 
wavenumber m, so that the normalized equal-time density-density correlations follow 



im) {q^) v^ qI 
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